%load_ext autoreload
%autoreload 2

import warnings
from functools import partial

from optimization import jg_opt
from pisa.administrative_area import AdministrativeArea
from pisa.facilities import Facilities
from pisa.population import WorldpopPopulation
from pisa.population_served_by_isopolygons import get_population_served_by_isopolygons
from pisa.visualisation import (
    plot_facilities,
    plot_isochrones,
    plot_population,
    plot_population_heatmap,
)

INFO:matplotlib.font_manager:Failed to extract font properties from /usr/share/fonts/truetype/noto/NotoColorEmoji.ttf: Can not load face (unknown file format; error code 0x2)

INFO:matplotlib.font_manager:generated new fontManager

import os

from dotenv import load_dotenv

# load environment variables from an `.env` file in the local root directory
load_dotenv()

CBC_SOLVER_PATH = os.getenv(
    "CBC_SOLVER_PATH"
)  # path to the cbc executable (e.g. /opt/homebrew/bin/cbc)

Pisa Showcase, Example Notebook

Define Administrative Area

Let op: the administrative area is a country because of the package that we use for fetching the data (gadm).

The naming is confusing: we want the administrative area to be a subset of a country (in this case, the administrative area should be Baucau).

Should be easy to fix once we change gadm to pygadm (see issue #59 on GitHub)

timor_leste = AdministrativeArea(country_name="Timor-Leste", admin_level=1)

# these are the boundaries of Baucau
# type: Polygon
baucau = timor_leste.get_admin_area_boundaries("Baucau")

INFO:pisa.administrative_area:Validating country name: Timor-Leste

INFO:pisa.administrative_area:Country name 'Timor-Leste' validated successfully

INFO:pisa.administrative_area:Retrieving boundaries of all administrative areas of level 1 for country Timor-Leste

  0%|          | 0.00/6.25M [00:00<?, ?B/s]

  1%|          | 32.8k/6.25M [00:00<00:27, 227kB/s]

  3%|▎         | 209k/6.25M [00:00<00:07, 784kB/s] 

 14%|█▍        | 889k/6.25M [00:00<00:02, 2.46MB/s]

 45%|████▌     | 2.83M/6.25M [00:00<00:00, 6.71MB/s]

100%|██████████| 6.25M/6.25M [00:00<00:00, 9.14MB/s]

Facilities

Get existing facilities (in our case, hospitals) from OSM

# Suppress user warning about geometry in geographic CRS. Centroid is calculated over a single facility (e.g. a hospital),
# so distances are very small and projection isn't necessary
warnings.filterwarnings(
    "ignore",
    message="Geometry is in a geographic CRS. Results from 'centroid' are likely incorrect",
    category=UserWarning,
)

hospitals_df = Facilities(admin_area_boundaries=baucau).get_existing_facilities()

INFO:pisa.facilities:Retrieving existing facilities with tags {'amenity': 'hospital', 'building': 'hospital'} using OSM.

INFO:pisa.facilities:Successfully retrieved existing facilities from OSM.

plot_facilities(hospitals_df, baucau)

Make this Notebook Trusted to load map: File -> Trust Notebook

Estimate potential locations for new facilities

potential_hospitals_df = Facilities(
    admin_area_boundaries=baucau
).estimate_potential_facilities(spacing=0.05)

Get population

In this example, from WorldPop

# TODO: why "ID" column? extra index for what reason? Maybe becomes apparent later?

population_gdf = WorldpopPopulation(
    admin_area_boundaries=baucau, iso3_country_code=timor_leste.get_iso3_country_code()
).get_population_gdf()

population_gdf.head()

	ID	longitude	latitude	population	geometry
0	0	126.14917	-8.67333	0.10	POINT (126.14917 -8.67333)
1	1	126.14917	-8.67250	0.09	POINT (126.14917 -8.6725)
2	2	126.15000	-8.67333	0.10	POINT (126.15 -8.67333)
3	3	126.15000	-8.67250	0.09	POINT (126.15 -8.6725)
4	4	126.15000	-8.67167	0.09	POINT (126.15 -8.67167)

plot_population_heatmap(population_gdf, baucau)

Make this Notebook Trusted to load map: File -> Trust Notebook

plot_population(population_gdf, baucau, random_sample_n=1000)

Make this Notebook Trusted to load map: File -> Trust Notebook

Calculate isopolygons

Here we make some choices:

• distance type

• distance values

• mode of transport

Valid values for constants are in the script pisa.constants

DISTANCE_TYPE = "length"

DISTANCE_VALUES = [2000, 5000, 10000]

MODE_OF_TRANSPORT = "driving"

Using OSM

For this, we need to get the road network from osmnx

from pisa.osm_road_network import OsmRoadNetwork

road_network = OsmRoadNetwork(
    admin_area_boundaries=baucau,
    mode_of_transport=MODE_OF_TRANSPORT,
    distance_type=DISTANCE_TYPE,
).get_osm_road_network()

INFO:pisa.osm_road_network:OSM road network set with parameters network_type 'drive' and distance_type 'length'

Calculate isopolygons for existing facilities

from pisa.isopolygons import OsmIsopolygonCalculator

isopolygons_existing_facilities = OsmIsopolygonCalculator(
    facilities_df=hospitals_df,
    distance_type=DISTANCE_TYPE,
    distance_values=DISTANCE_VALUES,
    road_network=road_network,
).calculate_isopolygons()

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

isopolygons_existing_facilities.head()

	ID_2000	ID_5000	ID_10000
osmid
1777614876	POLYGON ((126.59878421714033 -8.54764218472667...	POLYGON ((126.59761276961231 -8.57116537023302...	POLYGON ((126.56068305448515 -8.56751807760911...
1777614896	POLYGON ((126.65907104033573 -8.62868878467725...	POLYGON ((126.6311023795325 -8.646327183432366...	POLYGON ((126.63106042126435 -8.64641589673682...
1777615024	POLYGON ((126.59499222126435 -8.47806089673682...	POLYGON ((126.55950911149382 -8.48560952854904...	POLYGON ((126.57412141532275 -8.51524054033573...
1777615035	POLYGON ((126.68354761081066 -8.63974323564052...	POLYGON ((126.68357328721343 -8.63970069138097...	POLYGON ((126.6311023795325 -8.646327183432366...
1777615036	POLYGON ((126.6803646 -8.578027, 126.680359784...	POLYGON ((126.6650524795325 -8.602129783432366...	POLYGON ((126.65894556961231 -8.62895407023302...

plot_isochrones(isopolygons_existing_facilities, baucau)

Make this Notebook Trusted to load map: File -> Trust Notebook

Calculate isopolygons for potential facilities

isopolygons_potential_facilities = OsmIsopolygonCalculator(
    facilities_df=potential_hospitals_df,
    distance_type=DISTANCE_TYPE,
    distance_values=DISTANCE_VALUES,
    road_network=road_network,
).calculate_isopolygons()

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry

/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:402: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  disks = nodes_gdf.buffer(node_buffer).geometry
/home/runner/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/pisa/isopolygons.py:407: UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect. Use 'GeoSeries.to_crs()' to re-project geometries to a projected CRS before this operation.

  cylinders = edges_gdf.buffer(edge_buffer).geometry

Prepare optimization data

Preparing the variables that go into the cell below

population_served_current = get_population_served_by_isopolygons(
    grouped_population=population_gdf, isopolygons=isopolygons_existing_facilities
)

current = {DISTANCE_TYPE: population_served_current}

population_served_potential = get_population_served_by_isopolygons(
    grouped_population=population_gdf, isopolygons=isopolygons_potential_facilities
)

potential = {DISTANCE_TYPE: population_served_potential}

population_count = population_gdf.population.values

BUDGET = [
    5,
    20,
    50,
]  # budget for the optimization in terms of how many locations can be built

cbc_optimize = partial(jg_opt.OpenOptimize, solver_path=CBC_SOLVER_PATH)

values, solutions = jg_opt.Solve(
    household=population_count,
    current=current,
    potential=potential,
    accessibility=DISTANCE_TYPE,
    budgets=BUDGET,
    optimize=cbc_optimize,
    type="ID",
)

WARNING:pyomo.opt:Failed to create solver with name '_cbc_shell':
Failed to set executable for solver cbc. File with name=D:\EiriniK\Downloads\amplbundle.mswin64\ampl.mswin64\cbc.exe either does not exist or it is not executable. To skip this validation, call set_executable with validate=False.
Traceback (most recent call last):
  File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/base/solvers.py", line 148, in __call__
    opt = self._cls[_name](**kwds)
  File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/solvers/plugins/solvers/CBCplugin.py", line 88, in __init__
    super(CBCSHELL, self).__init__(**kwds)
  File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/solver/shellcmd.py", line 66, in __init__
    self.set_executable(name=executable, validate=validate)
  File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/solver/shellcmd.py", line 115, in set_executable
    raise ValueError(
ValueError: Failed to set executable for solver cbc. File with name=D:\EiriniK\Downloads\amplbundle.mswin64\ampl.mswin64\cbc.exe either does not exist or it is not executable. To skip this validation, call set_executable with validate=False.

WARNING:pyomo.opt:Failed to create solver with name 'cbc':
Attempting to use an unavailable solver.

The SolverFactory was unable to create the solver "_cbc_shell"
and returned an UnknownSolver object.  This error is raised at the point
where the UnknownSolver object was used as if it were valid (by calling
method "set_problem_format").

The original solver was created with the following parameters:
	executable: D:\EiriniK\Downloads\amplbundle.mswin64\ampl.mswin64\cbc.exe
	type: _cbc_shell
	_args: ()
	options: {}
Traceback (most recent call last):
  File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/base/solvers.py", line 148, in __call__
    opt = self._cls[_name](**kwds)
  File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/solvers/plugins/solvers/CBCplugin.py", line 50, in __new__
    opt.set_problem_format(ProblemFormat.cpxlp)
  File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/base/solvers.py", line 113, in __getattr__
    self._solver_error(attr)
  File "/home/runner/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/base/solvers.py", line 116, in _solver_error
    raise RuntimeError(
RuntimeError: Attempting to use an unavailable solver.

The SolverFactory was unable to create the solver "_cbc_shell"
and returned an UnknownSolver object.  This error is raised at the point
where the UnknownSolver object was used as if it were valid (by calling
method "set_problem_format").

The original solver was created with the following parameters:
	executable: D:\EiriniK\Downloads\amplbundle.mswin64\ampl.mswin64\cbc.exe
	type: _cbc_shell
	_args: ()
	options: {}

---------------------------------------------------------------------------
RuntimeError                              Traceback (most recent call last)
Cell In[20], line 9
      1 BUDGET = [
      2     5,
      3     20,
      4     50,
      5 ]  # budget for the optimization in terms of how many locations can be built
      7 cbc_optimize = partial(jg_opt.OpenOptimize, solver_path=CBC_SOLVER_PATH)
----> 9 values, solutions = jg_opt.Solve(
     10     household=population_count,
     11     current=current,
     12     potential=potential,
     13     accessibility=DISTANCE_TYPE,
     14     budgets=BUDGET,
     15     optimize=cbc_optimize,
     16     type="ID",
     17 )

File ~/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/optimization/jg_opt.py:362, in Solve(household, current, potential, accessibility, budgets, optimize, type)
    360 I = np.unique(list(IJ.keys()))
    361 J = np.unique(np.concatenate(list(IJ.values())))
--> 362 optimization = optimize(
    363     household,
    364     I,
    365     J,
    366     IJ,
    367     [max(budgets)],
    368     parsimonious=True,
    369     maxTimeInSeconds=5,
    370     mipGap=1e-15,
    371 )
    372 optimization["nof"] = [len(s) for s in optimization.solution]
    373 coverage = (optimization.value / household.sum() + percent_covered).to_frame()

File ~/work/Public-Infrastructure-Service-Access/Public-Infrastructure-Service-Access/optimization/jg_opt.py:126, in OpenOptimize(w, I, J, IJ, budget_list, parsimonious, maxTimeInSeconds, mipGap, trace, solver, solver_path)
    124 result.at[p, "modeling"] = pc() - start
    125 start = pc()
--> 126 solver_result = solver.solve(M, tee=trace)
    127 result.at[p, "solving"] = pc() - start
    128 result.at[p, "value"] = int(
    129     np.ceil(M.weighted_coverage() - np.finfo(np.float16).eps)
    130 )

File ~/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/base/solvers.py:99, in UnknownSolver.solve(self, *args, **kwds)
     97 def solve(self, *args, **kwds):
     98     """Perform optimization and return an SolverResults object."""
---> 99     self._solver_error('solve')

File ~/.cache/pypoetry/virtualenvs/gpbp-OIyr1Vno-py3.10/lib/python3.10/site-packages/pyomo/opt/base/solvers.py:116, in UnknownSolver._solver_error(self, method_name)
    115     def _solver_error(self, method_name):
--> 116         raise RuntimeError(
    117             """Attempting to use an unavailable solver.
    118 
    119 The SolverFactory was unable to create the solver "%s"
    120 and returned an UnknownSolver object.  This error is raised at the point
    121 where the UnknownSolver object was used as if it were valid (by calling
    122 method "%s").
    123 
    124 The original solver was created with the following parameters:
    125 \t"""
    126             % (self.type, method_name)
    127             + "\n\t".join("%s: %s" % i for i in sorted(self._kwds.items()))
    128             + "\n\t_args: %s" % (self._args,)
    129             + "\n\toptions: %s" % (self.options,)
    130         )

RuntimeError: Attempting to use an unavailable solver.

The SolverFactory was unable to create the solver "cbc"
and returned an UnknownSolver object.  This error is raised at the point
where the UnknownSolver object was used as if it were valid (by calling
method "solve").

The original solver was created with the following parameters:
	executable: D:\EiriniK\Downloads\amplbundle.mswin64\ampl.mswin64\cbc.exe
	type: cbc
	_args: ()
	options: {'threads': 8, 'sec': 5.0, 'allowableGap': 1e-15}

values

	10000	5000	2000
5	0.385968	0.256445	0.126371
20	0.393635	0.266618	0.133419
50	0.393635	0.268916	0.134770

solutions

	10000	5000	2000
5	[27, 43, 16, 10, 17]	[27, 22, 25, 9, 43]	[39, 43, 25, 7, 22]
20	[27, 43, 16, 10, 17, 22, 18, 35, 37, 14, 38, 0...	[27, 22, 25, 9, 43, 16, 18, 42, 50, 38, 39, 28...	[39, 43, 25, 7, 22, 27, 24, 26, 13, 32, 10, 16...
50	[27, 43, 16, 10, 17, 22, 18, 35, 37, 14, 38, 0...	[27, 22, 25, 9, 43, 16, 18, 42, 50, 38, 39, 28...	[39, 43, 25, 7, 22, 27, 24, 26, 13, 32, 10, 16...